Entanglement Entropy for Disjoint Subsystems in 
XX Spin Chain 



B.-Q. Jin* and V. E. Korepin^ 

* College of Physics and Electronic Information Engineering, Wenzhou University, 
Wcnzhou 325035, P.R. China 

C.N. Yang Institute for Theoretical Physics, State University of New York at Stony 
Brook, Stony Brook, NY 11794-3840, USA 

E-mail: j inbqOwzu . edu . on , korepinQinst i . phy s i cs . sunysb . edu 

Abstract. Fishcr-Hartwig formula has been successful applied to describe the von 
Neumann and Renyi entropies of a block of spins in the ground state of XX spin chain. 
It was based on a determinant representation. In this paper, we generalize the free 
fermion method to obtain an exact formulation for the entropy of any finite subsystem 
in XX spin chain. Based on this, we derive a determinant representation of the entropy 
of multiple disjoint intervals in the ground state of XX model. 



PACS numbers: 02.30.Ik, 03.65.Ud, 05.30.Rt, 05.50.+q, 64.60.De 



Submitted to: J. Phys. A: Math. Gen. 



Entanglement Entropy for Disjoint Subsystems in XX Spin Chain 



2 



1. Introduction 



There is great interest to quantify the entanglement in extended systems in recent 
years since its abihty to detect the scahng behavior in proximity of quantum critical 
points. Both von Neumann entropy and Renyi entropy play a key role in these 
developments and are so-called entanglement entropies as the measures of entanglement. 
In the earher studies of quantum spin chain models, one usually divides the system 
into two subsystems with each only containing contiguous lattice sites and finds the 
entanglement entropies of any subsystem when the whole system is in a pure quantum 
state [H 121 El in El E]- Recently, there is also growing interest to find the entropy 
of a non-contiguous subsystem motivated by the search for phase-transition indicators 
[3 El El [im im 1121 US] and its relation to conformal field theory [3 El E]- 

The physical system we consider is XX quantum spin chain with lattice index 
{1, 2, ■ ■ ■}. The Hamiltonian for XX quantum spin chain can be written as 



(1) 



ra=l 



All these lattice sites are separated into two sets denoted by Ca and Cb respectively 
and we consider the entropy of subsystem A of spins located on set Ca when the spin 
chain in the ground state. 

Following [H], we define two sets of Majorana operators 



C21-1 = I n < and C2i= n < I 

in<i / \n<l 



on each site of the spin chain and 



(2) 



C21-1 =1 n < and C2/ = n <\ (3) 

on each site of set Ca- With the help of the first set of Majorana operators defined by 
(|2| , XY quantum spin chain has been exactly solved [HI [15] and following correlations 
have been obtained 

{GS\c,n\GS) = 0, {GS\cmCn\GS) = + i (B^.)^,. (4) 

Here \GS) denotes the ground state of system and matrix Bx, can be written in a block 
form as 

/Ho n_i . . . u^_L \ 
Hi Ho : 



B, 



(5) 



/ 



with 



n. 







1 / 


V -9-1 




' ''=2^ Jo 



(6) 
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For XX model, 

f 1 for -kp < 9 < kp 

9(9) = I 7 

(-1 foi kF<9< {271 -kp) 

and kp = arccos(|/i|/2). Therefore, we have 

{GS\c2m-lC2n\GS) = - {G S\c2nC2m-l\G S) = I Qm-n = ''i- Qn-m- (8) 

To find the entropy of spins located on set Ca, the essential step is to find the 
reduced density matrix for these spins. From |2j, one can express this reduced density 
matrix with the second set of Majorana operators defined by (|3| and their multiphcation 
terms. Similarly as suggested by [1], if one can find a set of fermion operators 
which is linearly equivalent to the set of {cj} (i.e. hi and 6^ are linear combination of 
Majorana operators {cj} and transform matrix is of full rank) to satisfy 

{GS\hi\GS) = {GS\b+\GS) = 0, (9) 
{GS\bibj\GS) = {GS\btb+\GS) = (10) 

and 

{GSlkbpGS) = {GS\kbt\GS)6,,„ (11) 
then one can find the reduced density matrix 



PA = n [{GS\bth\GS)bth + {GS\hbt\GS)hbt). (12) 

The reduced density matrix for some specific subsystems of XY model has also been 
discussed in [161 El [H] • 

In the next section, we will first try to find this set of fermion operators {bi,bf} 
satisfying ([9]), (10) and (11), then obtain the reduced density matrix and its eigenvalues, 
finally get the von Neumann entropy and Renyi entropy for any finite subsystem. 



2. Exact formulation for entropy of any finite subsystem 



In order to find (or to find whether there exists) a set of fermion operators {bi,b^} 
satisfying (|9|, (|lO|) and (|lT|) as done in [H El |3], one has to find all correlations of 
{GS\cm Cn\GS). In order to find these correlations, we first express Majorana operators 
{cj} in terms of Majorana operators {cj}. From definitions ^ and ([s]), we have 

a' 



"i C2n-lC2n ^ud 
C21-1 = 



[ n 0C2«-1, 



C2l = i n Oc2i- 

n<l,n&CB 



Therefore, 



C2/-1 = n (-iC2n-lC2n)c2/-l and C2l = ]J ("1 C2n-lC2n)c2Z. 



(13) 



(14) 
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With the help of Wick Theorem and Q , we have 

{GS\CmCn\GS) = {GS\Cm Yl ["1 C2fc-lC2fc] C„ | GS") , (15) 

fce£sn£[(m+i)/2][(„+i)/2] 

where Cmn denotes the set of lattice sites between (but not including) m-th and n- 
th sites, and [x] means the integer part of x. Therefore, set C[(^m+i)/2][{n+i)/2] includes 
all lattice sites between (but not including) lattice sites on which Majorana operators 
Cm and c„ are defined. The right-hand side of the above identity for correlations 
{GS\cm Cn\GS) can be further expressed by correlations {GS\cpCg\GS) according to Wick 
Theorem. Let's assume there is K lattice sites inside of the set Cb H 'C[(^m+i)/2][(n+i)/2] 
with lattice index {Di, D2, • • • , D^}, then by brute force one can find 

{GS\c2p-i52,\GS) = idet{Tp,) (16) 

and 

{GS\c2p~C2,\GS) = {GS\c2p-i~C2,-i\GS) = S„. (17) 

Here {K + 1) x [K + 1) matrix Tpq is Toeplitz-like defined with row index i G 
{p, Di, D2, ■ ■ ■ , Dk}, column index j G {q, Di, D2, • ■ ■ , D^}, and matrix element in z-th 
row and j-th column by 

iTp,h = 9^-, = ^ r e-'^^-^^'gi9) (18) 



-PUzj 27r Jo 

and g{9) in ([7|. Interchanging row and column of Tpg and noticing gi = g_i for XX spin 
chain, one can find that 

det(Tp,) = det(T,p). (19) 

Now let us define 

ai = ^{c2i-i - ic2i) then = ^ {C21-1 + i C21) if / G £a- (20) 



From (35), we have 

{GS\di\GS) = {GS\dl\GS) = 0. (21) 
Let us define matrix A [18j through 

{GS\dm d+\GS) = ^ + i^^. (22) 



Then from (17) we can find that 



i^)mn = ^ ((^^152^-1 ~C2n\GS) + {GS\c2n-l C2m\GS)) . (23) 



From (16) and (19), we obtain that 

i^)mn = -\ [det(T^„) + det(r„^)] = - det(T„„). (24) 
Hence, in XX quantum spin chain, A is a real symmetric matrix with 

(A)_ = -det(T„„). (25) 
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Here Tmn is Toeplitz-like matrix defined in (|18|). Similarly, we have 
{GS\aiam\GS) = {GS 



{GS\dld+\GS) = {GS 



2 2 

C2I-I + iC2l C2m-1 + iC2m 



GS) = 0. 



GS) = 0. 



(26) 



(27) 



2 2 

Now let us take an example, in which set Ca contains the lattice site 1 and lattice site 
3, and work out its matrix A explicitly. For this case, A is a two by two matrix with 
row (and column) index taking {1,3}, i.e. 

V ^31 ^33 / 



(2J 



All matrix elements of above matrix can be obtained by (25). We have 



A 



33 



-go- 



(29) 



We also know that matrix element Ai^ is related to the matrix T13, which is defined in 



13; 



(18) with row index taking {1,2} and column index taking {3,2}. Hence, we have 

T,s = ( = (30) 

V 9-1 9o J \ 9i 9o J 

Therefore, we obtain 

Ai3 = A^i = gf - g2 go- (31) 

Since matrix A is a real symmetric matrix, it's diagonalizable with an orthogonal 
transformation, i.e. there exists an orthognal matrix V satisfying 

VAV'^ = Arf, (32) 



where A^ is diagonal matrix. Let's denote (V^)„„ = fmn 
bi = J^mVimdm, then we have bf = Em Vimd:^- Therefore, 

{GS\bib'j'\GS) = J2vim{GS\amd^\GS)vj„ = 



(A. 



1 + i^i . 
2 ^ 



i^m^mn and define 



Similarly, from (26) and ( p7| ), we have 

{GS\bibj\GS) = {GS\btb+\GS) = 0. 

With the help of Wick Theorem and (|4]) , we also have 
{GS\c^\GS) =0. 

Therefore, 

{GS\bi\GS) = {GS\b+\GS) = 0. 



(33) 

(34) 
(35) 
(36) 



Thus we find a set of fermion operators bi and b^ satisfying (9), (10) and (11) through 
the diagonalization of matrix A and obtain 

l + l^i r_L r I - l^i 



{GS\b, bt\GS) 



{GS\btk\GS) 



(37) 
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where {ui} are the eigenvalues of matrix A. Once we get this set of {i^i}, the von 
Neumann entropy of subsystem A can be expressed as 



siPA) = E 

i 

and Renyi entropy 

SaiPA) = E f 



m 



In 



a 



2 



— m 



1 - I'i 



(38) 



a ^ 1 and a > 0. (39) 



The way presented above can be used to exactly albeit numerically determine the 
entropy of system A containing several disjoint intervals. 



3. Determinant representation for entropy of any subsystem 



In order to develop a way of analytical treatment as in [2] , let us define 

A = A/-A, D(A) = detA, (40) 

where / is the identity matrix. Then we will have the determinant representation for 
the von Neumann entropy of subsystem A 

1 <p dAe(l + e,A) ^ 



S{pa) = 
with 

e(x, u) = 
and Renyi entropy 

Sa{pA) -- 

with 



lim — / dAe(l + e,A)4^1nD(A) 
e^o+ 27ri Jv dA 



X + V X + V X — V X — P 
In In 



(41) 
(42) 



Min^i dAe„(l + e,A)-^lnZ}(A) 
e^o+ 2tci Jf' dA 



1 — a 



In 



X + z/\ " ( X — V 



(43) 
(44) 



Here the contour F' is depicted in Figure [T| which encircles all zeros of -D(A). 

For the case of set La only containing contiguous sites P, A can be expressed as 
Toeplitz matrix with Toeplitz operators 

QA{Q) = \^g{Q) (45) 

where g{Q^ is defined in ([t]). With the help of Weiner-Hopf factorization and Fisher- 



Hartwig conjecture, integeration of (41) and (43) has been explicitly worked out for this 
case [2J. 

Although it's harder to give an explicit expression for A in the case of set La 
containing non-contiguous sites, A is of the form 



/ A. 



A 



"2 



(46) 



A. 
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Figure 1. Contours V (smaller one) and T (larger one). Bold lines (— oo,— 1 — e) 
and (1 + e, oo) are the cuts of the integrand e(l + e, A). Zeros of D]^{X) are located on 
the bold line (—1, 1). The arrows indicate the directions of integrations, and r and R 
are the radius of the circles. % 



with row (and column) index taking 1, 2, ■ ■ ■ , rii — 1, rii, ni + fci + 1, rii + /ci + 2, ■ ■ ■ , ni + 
ki + n2,ni + ki + n2 + k2 + 1, ■ ■ ■ , if the first ni lattice sites belong to Ca, then the 
next ki lattice sites belong to £b, then the next n2 lattice sites belong to Ca, and so 
on. These n, x sub-matrix A„. are Toeplitz matrix with generator ^a(^) defined in 
(45), and the elements of off-diagonal block (A) = — (A)., related to the determinant 



oTToeplitz-like matrix as seen in (25). So far, we got the determinant representation 
for the subsystem A containing several disjoint intervals. To have analytical treatment 
further, it will be necessary to tackle with the determinant of Toeplitz-like matrix T^n 
defined in (18). 



4. An application to two disjoint interval subsystem 

Now let us come to a simple but non-trivial case [7], in which subsystem A contains two 
disjoint intervals with with equal length m and distance m, i.e. lattice sites with index 
{1, 2, ■ ■ • , m} (denoted as Ai) and {2m -|- 1, 2m + 2, - ■ ■ , 3m} (denoted as A2) belong 
to subsystem A, but lattice sites with index {m -|- 1, m -|- 2, ■ ■ ■ , 2m} do not belong to 
subsystem A. We also take the external field h = in the model Hamiltonian ([T]). 
Therefore, the ground state of this system is in critical phase. We want to find the 
mutual information [71 18] between these two intervals, which is defined as 

Iai:A2 = Sai + Sa2 — Sa, (47) 
when the system is in the ground state. 
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From (25), we have the matrix A for subsystem A expressed as a block matrix 

^ A2,l A2,2 / ' 



A 



(4J 



Here related blocks can be expressed as 



Ai,i = A 



2,2 



9o -9-1 ••• -5-1- 



Ai,2 = A^i 



~9i 

\ —9m-l 
( All Al2 
A2I A22 



-9o 



(49) 



-90 J 



\ A 



ml 



A^ \ 



Amm J 



(50) 



with 



A. 



9i—j—2m 9i—'m—l 9i—'m—2 ■ ■ ■ 9i—'2m 

91- j-m 90 9-1 ■ ■ ■ 9l-m 

92- j-m 9l 9o '■ 



and from (|6| 



9-j 



-IV 



9j 



sm 



9m-l 



JJK 

2 



9o 



(51) 



(52) 



Finding all eigenvalues of matrix defined in (48) and substituting them into (38) 



one can find Sa- Similarly, one can obtain Sa^ = Sa2 by finding all eigenvalues 



of matrix defined in (49) and substituting them into (38). The mutual information 



between these two intervals is presented in Figure [2] From Figure [2| we find that the 
mutual information between two intervals will vanishes possibly when both the length 
of intervals and the distance between two intervals goes into infinite in the same scale 
even the system is in critical phase. It coincides the one obtained in in the case that 
the ratio of interval length to chain length goes to zero. 



5. Summary 

We have generalized the numerical calculation method for the entropy of one interval, 
which first appeared in [I], to the case of multiple disjoint intervals in the ground state 
of XX model. We use this formulation to find the mutual information between two 
intervals with the separation and the length of interval same. We confirm that the 
mutual information will vanishes possibly when the length scale goes into infinite even 
the system is in critical phase. We also derive the determinant representation of the 
entropy of multiple disjoint intervals for further analytical treatment. 
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0.04 



0.03 



0.01 0.02 0.03 0.04 



Figure 2. Iji-^.yi^(y-axis) as function of 1/m (x-axis). All points except (0,0) 
are calculated by taking m = 21,40,41,63,80,81,160,161,189,320,321,567,640,641 
respectively. Here,lAi:A2 the mutual information between two intervals, m is ttie 
length of intervals and also the distance between two intervals. % 
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